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SAFETY, RELIABILITY AND STRUCTURAL DESIGN* 


SYNOPSIS 


the safety analysis any structure the expected operational life for which 
is, should designed, represents important design parameter, the 
significance which not sufficiently understood. The concept design for 
infinite life unrealistic that design for zero probability failure, 
since the design life and the probability failure are related through the 
liability function” the designed expressing the probability 
survival function operational life, this function permits rational ap- 
proach the safety analysis for the criteria ultimate 
load failure and fatigue, well forthe combination both criteria. 

reality the safety factor shown Statistical distribution function 
the quotient resistance and load considered statistical variables. 
engineering significance can, therefore, attached specific single value 
this factor other than “measure location” (mean, median, made) 
the statistical distribution, However, the probability risk failure 
survival unique parameter the distributicn function the safety factor, 
and should replace the conventional concept safety rather considered 
its only rational quantitative Its generalization the 
“reliability function” provides the basis for integrated safety analysis 


Note.—Discussion open until August.1, 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Structural Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 87, No. March, 1961. 

This research was supported part the United States Air Force under Contract 
AF-33(616)-6288, monitored the Materials Laboratory Wright Air Development Divi- 
sion. 

prof. Civ. Engrg., Columbia Univ., New York, 
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structures permitting quantitative differentiation between the safety non- 
redundant and redundant 


INTRODUCTION 


The increasing severity operating conditions modern engineering struc- 
tures and the increasing complexity their structural action makes increas- 
ingly difficult such structures for assumedly infinite, but fact only 
unspecified, operational life without adversely affecting either their perform- 
ance, their cost, both. Moreover, the growing use materials high 
strength but unimproved fatigue performance increases the fatigue sensitivity 
such structures and requires the methods safety analysis 
which the criteria ultimate-load failure and failure under repeated load 
(fatigue proper alternating plasticity) are considered jointly. Such methods 
must obviously based quantitative measure safety that retains its 
validity under both failure conditions. With the aidof such safety, 
joint effect the combination failure criteria can adequately expressed. 

This measure the probability survival failure 
(“mortality”) the structure the termination the operational life for 
which has been designed or, the operational life has not been specified, 
function the (variable) operational life. For conditions ultimate load 
failure well the operational life expressed terms 
the expected number load repetitions throughout the service life the 
Thus the “reliability function” L(N) its complement, the “mor- 
tality function” P(N) L(N), expressing, respectively, the probabilities 
survival failure, functions the momentary “age” the structure, 
provide effective quantitative measures safety that are valid for either fail- 
ure criterion well for combination both them. 

Such combination most expediently based the concept the “risk” 
failure the “failure rate” r(N), which probability failure 
the 1)th load repetition those structures structural parts the 
initial statistical “population” considered that have survived load repetitions, 
and have thus reached the “age” definition, the “risk func- 
tion” r(N) the probability failure the statistical population the Nth 
load cycle p(N), referred the surviving percentage the population the 
Nth cycle L(N). Because the relation between the “probability density” p(N) 
and the “cumulative probability” reliability function L(N) 


the “risk function” can stated 


Therefore, integration 


- 

bik 

x 
is 
+: 

i 


statistical theory the riskfunction known intensity function, and its 
inverse the Mills 

While the reliability function L(N) purely statistical measure dis- 
persion, the risk function r(N) open physical (or engineering) interpre- 
tation, Under conditions ultimate load failure the “risk” failure is, 
definition, independent the “age” the structure and therefore constant, 
while the nature fatigue failure implies “rate” failure that in- 
creases with “age.” total “risk” under the combination the two failure 
criteria simply the sum the “risks” associated with the individual cri- 
Denoting the (constant) probability risk ultimate load fail- 
ure under single load application and the (increasing) risk 


logN 


FIG, 1.—RELIABILITY FUNCTIONS (N) FOR ULTIMATE LOAD FAILURE, 
FATIGUE, AND COMBINATION BOTH 


fatigue failure the N-th load application, the combined risk failure3 be- 
comes 


and, therefore, the associated reliability function 


for which the form the reliability function for fatigue failure 


depends the form the risk function (N) and vice versa (Fig. 1). 


1958, (a) 20; (b) pp. 75-155; (c) pp. 272-302; (d) pp. 255-272. 

“Prediction Fatigue Failure,” Freudenthal, Journal Applied Physics, 
Vol. 31, December, 1960. 
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and are based the simplifying assumption that and 
are independent, which implies that the possibility gradual reduction the 
ultimate carrying capacity spreading fatigue damage well that fa- 
tigue damage produced load-intensities ultimate failure-load range 
are both neglected. Without this assumption, must introducedas func- 
tion (N) considering the effect fatigue-crack propagation the 
“static” carrying capacity. addition would have determined 
from fatigue tests under variable load amplitudes including not only the range 
operational loads the “limit-load” serviceability (that might occur 
most once few times during the service life the structure), but also 
the extremely rare load-amplitudes exceeding this limit, but without causing 
failure single application. Such fatigue tests are not only impractical be- 
cause their duration, but the effects loadamplitudes that are 
occur during the expected life are, fact, negligibly small. However, the 
effect spreading fatigue damage can only neglected long 
function (N) rather than constant value must considered. 


SAFETY ANALYSIS FOR ULTIMATE LOAD CRITERION 


The probability failure under single load application refers 
statistical “population” nominally identical structures (statistical) 
carrying capacity each which has been subjected one application 
(statistical) load This represents the ratio structures that are expected 
fail, the probability any those structures fail, under the single 
load the “safety-factor” structure within the considered 
“population” defined the ratio its resistance the applied load 


This ratio different for each structure the is, fact, 
statistical variable with probability density p(v) and cumulative function P(v) 


p(v) dv, which completely determined the probability densities 


the carrying capacity and the load population 

The term denotes the carrying capacity the structure part the 
structure determined the relevant material properties and the geometric 
characteristics associated with the failure criterion, such areas, section 
moduli, static spans, well the effect connections (stress- 
concentration factors). has the same dimension the load which rep- 
resents the entity possible load intensities that might encountered the 
structures, even their occurrence quite unlikely within the period ser- 
vice for which the structure Whereas always statistical 
variable, there are types structures for which may represented 
non-statistical maximum load Smax. This because true service load 
high frequency occurrence, such loads ware-houses, bulk-and 
fluid-containers, cranes, railroad-bridges, and on, For such structures 
the distribution p(v) 


identical with the distribution divided the constant Smax. 


identical with the probability that P(1), which the ordinate 
the cumulative function P(v) The relationship 


determines the probability failure directly terms the probability den- 
sities and p2(S) without reference specific value the “safety fac- 
tor.” 

The definition the safety factor terms statistical distribution 
function p(v) implies that any value reduces the 
safety analysis structure structural part the determination, the 
basis Eq. 7a, the function p(v) and the integration specified Eq. 
the distribution p(v) unimodal, values within certain “central” range 
(around the mode) are the most likely, whereas the likelihood values out- 
range decrease with increases their distance from the central range. 
Specific values acquire meaning only they represent “measure 
location” the central range, such the mean mode median the 
distribution p(v). conveniently selected “measure location” can con- 
sidered nominal “central” safety factor representing the ratio the 
central values The selection this ratio will depend onthe charac- 
ter and p2(S). Thus, these functions are normal, their 
“central” values are the means Rand the values are logarithmic-normal, 
the central values, being the means (log and (log are the medians 
and the probability level 0.5, since log log log log 
the functions are their central values are the 
1/e) and (1/e), respectively. 

Two types structures must considered: 


one between independent statistical variables. can, therefore, obtained 
relatively simple statistical 

Structures which depends through the effect the flexi- 
bility and inertia. The ratio is, therefore, one dependent 
variables andcannot derivedunless dependence considered 
the statistical analysis. 


Structures the first type are those for which stress-analysis can based 
the undeformed shape structures). The intensity the forces for 
structures the second type depends significantly the deformational re- 
sponse the structure which, turn, affected its carrying capacity 
(“flexible,” “aero-elastic,” and “hydro-elastic” structures). Only structures 
the first type will considered. 

The principal feature rational approach the safety analysis such 
structures the consideration the safety factor statistical variable 
probability density p(v) determined the probability densities carrying 
iceability,”4 the (statistical) safety factor and probability re- 
fer the same refers toa criterion actual failure 
(fracture, plastic collapse) the associated (statistical) safety factor 


“Synopsis First Progress Report, the Committee Factors Safety,” 
Julian, Proceedings, ASCE, Vol. 83, No. July, 1957. 
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Eqs. and are based the simplifying assumption that and 
are independent, which implies that the possibility gradual reduction the 
ultimate carrying capacity spreading fatigue damage well that fa- 
tigue damage produced load-intensities ultimate failure-load range 
are both neglected. Without this assumption, must introducedas func- 
tion (N) considering the effect fatigue-crack propagation the 


carrying capacity. addition would have determined 


from fatigue tests under variable load amplitudes including not only the range 
operational loads the “limit-load” serviceability (that might occur 
most once few times during the service life the structure), but also 
the extremely rare load-amplitudes exceeding this limit, but without causing 
failure single application. Such fatigue tests are not only impractical be- 
cause their duration, but the effects loadamplitudes most unlikely 
occur during the expected life are, negligibly However, the 
effect spreading fatigue damage can only neglected long 
function (N) rather than constant value must considered, 


SAFETY ANALYSIS FOR ULTIMATE LOAD CRITERION 


The probability failure under single load application refers 
statistical “population” nominally identical structures (statistical) 
carrying capacity each which has been subjected one application 
(statistical) load This represents the ratio structures that are expected 
fail, the probability any one those structures fail, under the single 
load the “safety-factor” structure within the considered 
defined the ratio its resistance the applied load 


This ratio different for each structure the “population,” is, fact, 
statistical variable with probability density p(v) and cumulative function 


p(v) dv, which completely determined the probability densities 


the carrying capacity and the load population 

The term denotes the carrying capacity the structure part the 
structure determined material properties and the geometric 
characteristics associated with the failure criterion, such areas, section 
moduli, static spans, well the effect connections (stress- 
concentration factors). has the same dimension the load which rep- 
resents the entity possible load intensities that might encountered the 
structures, even their occurrence quite unlikely within the period ser- 
vice for which the structure Whereas always statistical 
variable, there are types structures for which may represented 
non-statistical maximum load Smax. This because true service load 
high frequency occurrence, such loads ware-houses, bulk-and 
fluid-containers, cranes, railroad-bridges, and on, For such structures 
the distribution p(v) 


Smax 


identical with the distribution divided the constant Smax. 


identical with the probability that P(1), which the ordinate 
the cumulative function P(v) The relationship 


determines the probability failure directly terms the probability den- 
sities and p2(S) without reference specific value the “safety fac- 
tor.” 

The definition the safety factor terms statistical distribution 
function p(v) implies that any value possible, reduces the 
safety analysis structure structural part the determination, the 
basis Eq. 7a, the function p(v) and the integration specified Eq. 
the distribution p(v) unimodal, values within certain “central” range 
(around the mode) are the most likely, whereas the likelihood values out- 
range decrease with increases their distance from the central range, 
Specific values acquire meaning only they represent “measure 
location” the central range, such the mean mode median the 
distribution p(v). conveniently selected “measure location” can con- 
sidered nominal “central” safety factor representing the ratio the 
central values (Ro/So The selection this ratio will depend onthe charac- 
ter the functions and Thus, these functions are normal, their 
“central” values are the means Rand the values are logarithmic-normal, 
the central values, being the means (log and (log S), are the medians 
and the probability level 0.5, since log log and log log 
the functions are their central values are the “characteristic” 
1/e) and (1/e), respectively. 

Two types structures must considered: 


one between independent statistical variables. can, therefore, obtained 
relatively simple statistical procedures, 

Structures which depends through the effect the flexi- 
bility and inertia. The ratio is, therefore, one dependent 
variables andcannot derivedunless dependence considered 
the statistical analysis. 


Structures the first type are those for which stress-analysis can based 
the undeformed shape (‘rigid” structures). The intensity the forces for 
structures the second type depends significantly the deformational re- 
sponse the structure which, turn, affected its carrying capacity 
(“flexible,” “aero-elastic,” and “hydro-elastic” structures). Only structures 
the first type will considered. 

The principal feature rational approach the safety analysis such 
structures the consideration the safety factor statistical variable 
probability density p(v) determined the probability densities carrying 
iceability,”4 the (statistical) safety factor and probability re- 
fer the same criterion. refers toa criterion actual failure 
(fracture, plastic collapse) the associated (statistical) safety factor 


“Synopsis First Progress Report, the Committee Factors Safety,” 
Julian, Proceedings, ASCE, Vol. 83, No. July, 1957. 
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and probability are those such failure, Obviously, the probability den- 
sity the load usually the same for both criteria. Since the central 
locations the central safety factors and, therefore, 

Considering two independent, non-negative, continuous statistical variables 
and with probability densities and and (cumulative) proba- 
bility functions and the joint distribution the two variables 
p(R,S) the product 


Introducing the quotient (R/S) yields 
and 
The distribution 


the joint distribution the quotient and the variable inte- 
grating over between the probability density p(v) obtained 
the “marginal” distribution 


The probability function 


Changing the order integration yields 


Eq. represents the relation between P(v) and the known probability functions 
Rand The probability failure is, therefore, 


noted that this equation closely related the Mellin transforms® the 
functions and 


“Methods Safety Analysis Highway Bridges,” Sixth Congress Internatl. Assn. 
Bridge and Structural Engrs., Stockholm, 1960, 655. 

“Some Applications the Mellin Transform Statistics,” Epstein, Annual 
Mathematical Statistics, Vol. 19, 1948, 370. 
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SAFETY ANALYSIS 


which represent the (s-1)t® moments the distributions and or, 
statistical terminology, the expected values “expectations” and 
Hence, the Mellin transform 


the distribution p(v) the expectation 


(s-1) the function Inversionof Eq. produces the distribution p(v) 
directly. 

The Mellin transform G(s) the probability function P(v) obtained with 
the aid Eq. 14; 


which the Mellin transform the probability function and 
that the probability distribution Inversion G(s) produces 
the probability function 

particularly simple illustration the procedure deriving p(v) from 
and without the use Eqs. and 14, arises when and 
are both logarithmic-normal with median values and Sand standard 

log logR logS 


follows that log the difference two normal distributions, also nor- 
mally distributed, with mean 


and standard deviation 


Therefore 


which 


i 
and 
log 
a 


The probability failure thus obtained directly from tables error func- 
tions 


ineffective the same time increases, This frequently the case when 
structural materials higher strength that show increased scatter the 
strength properties are 

Fig. Eq. presented for various combinations the ratios 
and which and denote the standard deviations and re- 


FIG, FAILURE FUNCTION CENTRAL SAFETY 
FACTOR (LOGARITHMIC-NORMAL DISTRIBUTIONS) 


spectively. Fig. presents the results the numerical evaluation Eq. 
under the assumption that the extremal distribution smallest val- 
(“Weibull distribution”) 


whereas the extremal distribution largest (“Frechet dis- 
tribution” 


for various combination the ratios and which and are the 

“characteristic values” and and are scale parameters that are inverse 

functions and Table presents the values for these curves, 
Introducing Eqs. and into Eq. the following expression obtained: 


which 
Therefore, the probability failure 


This integral can only evaluated numerical integration, The results have 
been presented Fig. for various ratios and 
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FIG, FAILURE FUNCTION CENTRAL SAFETY 
FACTOR (EXTREME-VALUE DISTRIBUTIONS) 


For non-statistical load and Eq. 27b can expressed 
closed form 
The parameters and ag, related tothe coefficients varia- 


tion and Table These relationships permit the 
selection parameters represent observations and 


| | | | | = 

| | | | | | 


Whether logarithmic-normal extremal distributions are applied 
particular safety analysis will depend the level the control imposed 
the production the structural material and the construction process, 
well the extent variability the 

Logarithmic-normal distributions have been found fairly well represent 
material parameters obtained under the conditions stringent plant control 
when the level quality control only fair poor, construction sites. 
With respect loads, distributions are likely represent 


TABLE 1,—PROBABILITY FAILURE FUNCTION 
CENTRAL SAFETY FACTOR 


(a) logarithmic normal distributions 


(b) Fig. extreme-value distributions 


ooocooocooco 


the complete load population, while extremal distributions are obtained anew 
load populationis derived therefrom, containing only loan intensities exceeding 
relatively high This new population will have higher cen- 


“Safety and the Probability Structural Failure,” Freudenthal, Transac- 
tions, ASCE, Vol. 121, 1956, 1337. 


(1) (2) (3) 
0.00 
0.10 
0.00 
0.10 
0.00 
0.10 
0.20 
0.20 
0.20 
0.30 
0.30 
4 


tral location than that the parent population while its standard deviation 
will appreciably lower. 


ULTIMATE LOAD CRITERION REDUNDANT STRUCTURES 


When the considered structure has single load-pathand consists parts 
“resistance” each which subject the full load (Fig. 4), failure 
the structure will occur when the weakest the parts fails (non-redundant 
Denoting the probability failure the individual 
member, the probability survival all members 


k=n 
Fk) 


therefore, the probability failure the structure 
k=n 


k=n 


The survivorship function for load repetitions 


which identical with that for the single-member structure, provided the in- 
dividual probabilities fulfill the condition Eq. 29. When all members 


TABLE 


logarithmic-normal extremal 


0.084 0.123 


are equal probability failure this condition has the form 1/n Pr. 
The individual members n-member chain-structure must, therefore, 
designed with higher “central” safety factor than equivalent single- 
member structure, order achieve the lower probability failure 

statically determinate structure members equivalent structural 
action chain-structure, the failure any one its members pro- 
duces failure the structure whole. Thus, the probability survival 
expressed Eq. 28, the probability failure Eq. 29, and the reliability 
function Eq. 30. The only difference that different members are subject 
different stress-intensity fromthe same external load-spec- 
trum, that the variation the values for individual members are likely 
much widerthan simple chain-structure nominally identical, iden- 


4 


tically loaded members. Thus, the probability failure the members 
are the failure range, may small negligible the sum Eq. 
29. The probability failure therefore composed the probabilities 
the critical members only rather than all members, 

When the considered structure multiple-load-path structure mem- 
bers and load path (Fig. 5), which (m/n) integer, the probability 


FIG.4.—MODEL NON- FIG. 5.—MODEL REDUNDANT 
REDUNDANT (SINGLE- (MULTIPLE-LOAD-PATH) 
LOAD-PATH) STRUC STRUCTURE 
TURE 


failure each the component chain structures subject (1/m)th the 
total load initially given the expression 


m 
4 _n k= — 
a n n ( 


Failure the structure itself implies consecutive failure the component 
chain-structures course which the probabilities pp, members 
the unbroken chains are suddenly increased one after another the chains 
This increase due the increase the the total load 
carried bythe The increase reduces the “central” safety 
factor the chain members from their initial value 


n = | m iy 

n 

: 


and constant. The increase from the initial values for the 
members the unbroken chains for the n(1 1/m) mem- 
bers, after the first chain broken, and subsequently and 
can obtained from the relevant diagrams Logarithmic normal 
distributions Rand are obtained from 25, introducing the step-wise 
decreasing values associated with each chain failure, 


m 


The probability failure the structure the probability consecutive fail- 
ure all component chains 


The probability survival the m-load-path structure 


compared Eq. for the single-load-path structure the same total 
number members, 

load-path structure members suggests that the addition alternative 
load-path reduces the probability under single load application sig- 
nificantly provided least the failure the first load-path 
(chain) already brings the probability PF,1 Close unity, little gained 
the addition alternative load 

The reliability function obtained similar combination the relia- 
bility functions the individual chains the different stages chain-failure: 


/ 


n/m 


(N) exp (-N 


The mortality the structure obtained from the mortality functions 
the component chains the form the product 


| 
ie 
5 
pe 
4 
FS 
ae 
q 
ie 
| 


that the reliability function the entire structure 


This compares with Eq. with for the single load path structure, 
For instance, for three-load-path structure 


Assume distributions carrying capacity andof load with 


Fig. reducing from 4.5 3.0 failure one chain reduces 
10-4, while failure the second chainfurther reduces the cen- 


With these values Eq. can written the form 


Lo(N) 


The increase the long-range reliability the structure compared that 
the single chain given the first term relatively small. For the second 
term only extends over the range values 101 whereas for 
the third term for Beyond these ranges the last two terms 
not contribute. 

The principal effect, however, that the three-load-path structure assures 
that the failure one-load-path still leaves structure with sufficiently high 
probability short range The reliability function the structure 
after the failure the first chain then reduced 


Hence, the structure reduced strength retains measure safety compa- 
rable that the undamaged structure, provided its operating life suit- 
able reduced, The condition 10-7 10-4 furnishes the relation, 
250, between the operational lives the undamaged and the damaged 
structure the same level “reliability” L(N). structure, designed 
for operational life 106 load applications can, after failure one-load- 
path, still operated for 103 load repetitions without reducing its 
“reliability.” 

The foregoing analysis illustrates the principle so-called “fail-safe” de- 
which quite extensively applied the design aircraft structures 
eliminate the catastrophic consequences failure one-load-path structures. 
The value represents the ratio between design life and maximum length 
inspection periods for structures which failures are not immediately ob- 
vious but have looked for during inspections. 


“Fatigue Aircraft Structures,” Freudenthal, Academic Press, New York, 
Y., 1957, 308. 


The risk failure associated with the reliability function multiple- 
load-path structure obtained applying Eq. Eq. 


For long operational livesN that r(N) the risk 
failure independent the redundancy the structure. However, for short 
lives the failure risk first approximation 


cause 


The result clearly illustrates the superiority the multiple-load-path struc- 
ture. The risk function approaches zero the number load applications 
tends towards zero, while approaches the constant value the single- 
load-path-structure only when the number load applications longer 
Thus, for relatively small number load applications, the multiple- 
load-path structure has much smaller probability failure and, therefore, 
afar higher safety comparative For struc- 
tures subject periodic inspection within intervals associated sharp- 
reduced probabilities failure, the provisionof redundancies ensures prac- 
tically absolute safety. 

indeterminate structure operationally roughly equiv- 
alent load-path chain structure, which each chain represents 
alternative independent determinate system obtained from the indeterminate 
system removing different sets redundant members, This multiple- 
chain structure subject (m+1) times the actual load intensity, that, 
initially, each chain carries the full load the structure, redundants con- 
secutively fail, the load the remaining redundants increases and does 
their probability failure through consecutive decrease vg. Eqs. 
are applicable, with the modification that the values PFmk 
may vary widely within each alternative determinate 

the preceding analysis, has been assumed that failure one redundant 
completely eliminates its structural action. the collapse analysis m-fold 
determinate elastic-plastic structures, this not the case. different ap- 
proach the ultimate-load analysis such structures required, for their 
failure involves the simultaneous attainment the limiting plastic force 
moment members sections rather than their consecutive elimina- 
tion. When the loading proportional and the location these critical mem- 
bers sections can assumed remain fixed, the probability failure 
plastic collapse (ultimate carrying capacity) equal the probability si- 
multaneous occurrence (m+1) plastic hinges members. However, where- 
the carrying capacity each hinge (plastic moment) each member 
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(plastic force) independent statistical variable, the applied moments 
forces these hinges members are not independently variable, but related 
the statistical external load the equations defining the collapse mechan- 
These equations arising from the application the principle virtual 
work? are usually the form 


force-distributions under conditions collapse. The distribution and 
Ro... are essentially determined the distribution the yield-limit. The 


aid which the probability failure under single load application, 
obtainedat and the survivorship function L(N) exp 


SUMMARY 


The safety factor structure defined terms the frequency dis- 
tribution the quotient resistance and load, and its relation with the prob- 
ability structural failure developed, The related concepts “structural 
reliability” and failure” are introduced and examined basis 
safety analysis for finite operational life non-redundant and redundant 
structures, 

Methods safety analysis based logarithmic-normal and extremal dis- 
tributions resistance and load are developed and the effect redundancy 
the risk and reliability function are 


“Plastic Analysis Structures,” Hodge, McGraw Hill Book Co., Inc., New 
York, Y., 1959, 32. 
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MATRIX ANALYSIS STRUCTURES CURVED SPACE 


Frank Baron,! ASCE 


SYNOPSIS 


general solution given matrix form for the analysis member 
curved segmental space and continuous betweentwo supports. The mem- 
ber may loaded any direction. 

The solution also yields general slope deflection equation which involves 
stiffness matrix, The stiffness matrix contains all stiffness and carry-over 
factors for the ends the member, 


INTRODUCTION 


matrix formulation given for the analysis structure which lies 
space and which continuous between two supports (Fig. 1). The structure 
may have any shape space andany variation incross sectionalong its length, 
may subjected loads any direction, moments about any axis, and 
linear angular discontinuities any given station the 
assumed that Hooke’s Law applicable and that the rotations and displace- 
ments the structure and small relative the over-all dimensions the 
structure, The definition smallangle considered acceptable the suc- 
ceeding analyses, Furthermore, the effects shear deformation andof changes 
length element structure are considered negligible. desired, 
adjustments the succeeding analyses may made for these 


Note,—Discussion open until August 1961, extend the closing date one month, 
written request must filed with the Executive Secretary, This paper part 
the copyrighted Journal the Structural Division, Proceedings the American So- 
ciety Civil Engineers, Vol, 87, No. March, 1961, 

Prof, Civ, Engrg., Univ. California, Berkeley, Calif, 
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The special structure which plane included inthis for- 
mulation, 

Signs.—The sign convention illustrated Fig. Three orthogonal axes 
(x, and are selected for reference, but part the structure need lie 
the planes defined these axes. right-hand rule used for these and 
succeeding orthogonal (James Michalos, ASCE, and the writer 
have previously left-hand rule.) 2(a), the structure loaded 
with external forces, and external moments, which are represented 
The components and along the and axes are positive 
when they are against the directions these axes, 

The moments and forces section the structure are designated 


with subscripts which are explained when introduced. define positive mo- 
ments and forces different sections the structure, forward path 


FIG, 


observer defined. The forward path observer indicated Fig. 
curved arrowalong the When the observer moves inthe direction 
this arrow, section the back the observer aback face and the sec- 
tion the front front face. Forces and moments are positive front 
face when they are the directions the axes. They are positive back 
face when they are against the directions the axes, 

The same sign convention used for rotations, displacements, and angle- 
changes, Fig. 2(b), angle-change rotation treated force along 
the axisof rotation, displacement along anaxis istreated amoment about 
the axis displacement. Consequently, computations geometry involving 
small angles are interchangeable with those statics.2,3 The reciprocal 

“Laterally Loaded Plane Structures and Structures Curved Space,” Frank 
Baron and James Michalos, Transactions, ASCE, Vol, 117, 1952, pp. 279-311. 

“Successive Corrections Pattern Thought,” Frank Baron, Numerical Me- 


thods Analysis Engineering, Chapter Grinter, Ed., Macmillan Company, 
1949, 
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MATRIX ANALYSIS 


this interchange possible, that is, computations statics can interchanged 
with those geometry which involves the definition small angle. 

The structure shown 2(b) the conjugate that shown Fig. 2(a). 
The loads conjugate structure are the angle-changes, dé, the real 


FIG, 


structure, and the forces and moments cross section the conjugate (or 
analogous) structure are the same, respectively, the rotations and displace- 
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tion the structure are designed 


with subscripts which are explained when introduced, 
Position vectors are designated 


Subscripts are used define the direction position For example, 


PROPERTIES ELEMENT SPACE 
m= My 


each end, The principal axes the cross sections are orthogonal and are 
designated and The unit vectors along these axes are designated 
and and are written matrix notation follows: 


Giz 


Ag = do (6) 


The vector product denoted bythe symbol holds because the prin- 
cipal axes are orthogonal. 

The scalar components along the and axes are givenby the sca- 
lar product which the asterisk denotes that transposed, Conse- 
quently, the magnitudes the angle-changes along the and axes can 
written 


0 0 dag3 
which the diagonal elements the matrix are given 


| 
4 
ee (9) 
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The diagonal elements the matrix are scalar and are the magnitudes 
the angle-changes about the axes caused moments one unit 
about the same axes, respectively. vector form, the angle-change caused 


The angle-change can also written 


which the matrix the right side the equation, Each ele- 
ment the matrix angle-change about axis designated the first 


FIG, 


subscript. The angle-change caused moment one unit about axis 
designated the second subscript. For example, daxy angle-change 
bout the x-axis caused moment one unit about the y-axis. Maxwell’s 
reciprocal relationship holds, and 


and 


T 
. 
g 
| 
(14 


Consequently, symmetric matrix, that is, 
From and 12, the following relationship obtained: 


Consequently, each element the matrix defined, For example, the 
element given 


GENERAL STRUCTURE SPACE BETWEEN TWO SUPPORTS 


Consider the structure shown Fig. The structure lies space and 
continuous betweentwo supports, the same asthat dealt with 
The present formulation differs from the preceding using matrix notation 
which lends itself simplified programming for digital 

The analysis the structure shown Fig, must fulfill the usual require- 
ments statically indeterminate structure, The requirements are follows: 


and 
which 


Note that the antisymmetric vector matrix contains the components the 
vector sucha way that the product r°P results the vector product the 
vectors and (Fig. 4). 

Geometry.—For any closed circuit, 


and 


provided that the angle-changes are small, For these cases, the requirements 
geometry are identical with those statics when angle-change rota- 
tion treated force, and displacement treated moment, The 
preceding statement the basis for the interchange 

Properties statical and geometrical relationships 
must consistent with the assumed properties materials, Hooke’s Law 
assumed, 

The structure shown Fig. may analyzed procedure consisting 
two steps follows: 


(1) Assume any statically possible distribution forces and moments along 
the structure, Unless this happens thecorrect distribution forces and 
moments, errors geometry will exist. 


(19) 
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(2) Subtract forces and moments whichpreserves the re- 
quirements statics but corrects for the errors geometry. 


matrix form, the procedure consists 
O = 09 (24) 


which designates the final moments and forces acting onany section the 
structure, and and designate the moments and forces which are associ- 
ated with steps and respectively, The rotations and displacements any 
section the structure are given 


The procedure automatically satisfies the requirements statics, remains 
satisfy the requirements geometry. 

The loading conditions for steps and are shown Fig They are 
shown for segment the structure with imaginary rigid urm extending 
from point the origin the and z-axes end the structure, 


FIG, 


imaginary rigid arm also extends from point end These arms are in- 
troduced for convenience closing the circuit. For step statics yields 


antisymmetric vector matrix containing the and components vector 
p- 
The requirements geometry, given Eqs. and 23, yield 
and 


which the summations are performed forthe closedcircuit pabp. The terms 
with the subscript are referred the origin 


i 
Mer 
Z 
= 
i 
J 
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Now let 
and 
and observe that (rp da)* da* Therefore, 


and Eqs. and can rewritten 

A mpi 1p” fpi = eee (35) 
and 

Mpi Jp fpi = Mpo +++ 42 (36) 


respectively. solution for mpj and fpj can obtained means 
and follows: Rewrite the equations matrix form 


which and are the first, second, and third matrices Eq. 37, 


and also 


respectively. Eq. when premultiplied yields 
Spi Kp ener eee ee ee (39) 


The solution given Eq. complete and involves the inversion 
matrix. general, this imposes serious problems the inversion 
can performed any computer with subroutine capable inverting 
matrix. However, the solution can obtained simpler form which 
involves the inversion matrices instead, For that purpose, premul- 
tiply Eq. and subtract the result from Eq. 36, thus obtaining 


Let 


and 
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and rewrite Eq. 
Jp fpi Moo (43) 


which when premultiplied yields 


The solution for obtained similar way. Premultiply Eq. 
and subtract the from Eq. 35, thus obtaining 


Let 
and 
and rewrite Eq. 
which when premultiplied yields 
Mpj = Ap Poo 
Therefore, Eqs. and can written 
The solution completed recalling that 
-r? 


and 
OP 


The preceding analysis stated terms any arbitrarily selected 
and axes, For computational purposes, suggested that the position 
selected somewhere near the average location the da’s the structure. 
the sensitivity succeeding inversions, The formal for doing this 
explained the next section. 

Centroidal Axes.— Consider the same structure shown Fig. but with the 
imaginary arm now extending from end positionthat will referred 
acentroid, For the generalized structure, more than one centroid re- 
quired matrix notation, the moments and forces shown Fig. 6(a) 
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(9) 


if 
N 
< 
N 
> 
a 
apa 
> 
| 
oO 
= 
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are designated 
Myci 
Mzci 


and 


yci 


The forces are along axes which, for convenience, are called centroidal axes, 
The positions these axes with reference the and axes are given 


* 0 Ty = Xy ; Xz 


Statics the rigid arm between and yields 
Oni = Co Opi (56) 


which the elements the matrix are arranged follows: 
=O 
Tr = 0 -Zy Yz 


Note that the matrix product moment with respect 
axis through and that 


The correction moments and forces any cross section the structure are 


given 


which 


Fig. 6(b) rotations about the and axes and displacements along 
the centroidalaxes are shown, Each axisof rotation selected correspond 
with that correction moment and each axis displacement corresponds 
with that correction Note that the matrix product yields 
the contribution rotation displacements along the centroidal axes, 


- 
(55) 
j 
(57) 
(59) 
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Therefore, considerations geometry for the rigid arm yield 


When desired, the inverse given 


The requirements geometry for closed circuit 


and 
dam= (65) 


can transformed the centroidal axes simple matrix manipulation, From 
the preceding development, 


Kp Spi = Apo (66) 
Substituting 
= Co” Foie (67) 
and 

Apo = ee eevee (68) 


into Eq. 66, and premultiplying the resulting equation yields 
which can rewritten 


Consequently, the elements can obtained from those means 


The centroidal preferred axes can located obtaining the elements 
Eq. 57. From Eq. 71, obtain 


ad A Tro > ly” (72) 
which the general form 
A = AB + Cc oeneteveeeeeeeeeeeeee#ee (73) 


The elements are selected such that diagonal matrix when the 
diagonal elements are these stipulations, the solution 


= 
A = 6 = I 0 6 = k A (61) 
and 
| 
es 
} 


March, 


which the diagonal elements are given 
matrix notation, the subscripts designate the diagonal elements ma- 
trix, 


With reference the centroidal axes, the correction moments and forces 


(76) 


(77) 


(78) 


(79) 

and 

The correction moments and forces any section the structure are given 
Eq. 59. 

Other general equation that holds for other axes besides 

the centroidal axes. For example, consider and axes that through 


from point For this case, Eq. transformed into 


and the antisymmetric matrix, 


Eq. can used for transferring the computations referred 
those referred point Furthermore, the same equations canbe used 


which 

which 
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for transferring the computations areas, first moments areas, and mo- 
ments inertias parallel axes. 


STIFFNESS MATRICES 


Stiffness and carry-over factors foruse amoment distributionor 
deflection method can obtained for the general member The fac- 


FIG, 


tors are assembled matrix form yield the fundamental slope-deflection 
relationships for the member, 
The member shown Fig. and subjected slopes and deflections, 
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(86) 


ends and respectively. loads are the member between ends and 
The moments and forces ends and are designated 


Mab 
fab 


which the first subscript each term designates the end the member 
which the term referred, Statics yields 


= I -Tab Ge 


which antisymmetric vector matrix containing the and com- 
ponents vector Statics the member between and yields 


The requirements geometry for the member between and yield 

and 

which can rewritten 

and 


these equations, and are the same kind those given Eqs. 
Eq. and are combined into the single matrix equation 


which 
“eee (96) 


and 
and 
. 
a 
and 
| 


The solution for given 


which when substituted into Eq. yields 


Assembling Eqs. and into single matrix equation yields 
which 
(101) 


(102) 


(103) 


Eq. 100 the fundamental slope-deflection equation for general member 
lying space, The equation relates the forces the ends the 
member the slopes and deflections the ends, means the properties 
the member. The properties are given which called the stiffness 
matrix the member, The stiffness matrix and involves the in- 

The solution for can obtained simpler form which involves the 
inversion matrices instead, The latter form also clarifies the phys- 
ical significance each sub-matrix and element For these purposes, 
rewrite Eqs. and 


which, for example, designates the moments and forces caused 
unit rotations and unit displacements when end fixed, Eq. 104 
simply statement proportionality between causes and effects for the given 
member, 
the same way the previous section, Eqs. and can operated 


a 

q 

atte 
Ree 
and 
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which and are the same kind terms those given Eqs. and 
46, From 89, 104, and 106, obtain 


[Ep)a Ep)b| = [Ka Ea)a = Ka Ea)a Kp ] (107) ; 
which completes the solution for the inverses matrices 


Values of m,, Volues of m,, 


FIXED END 
END 


MEMBER 


MEMBER 

108 


Stiffness and Carry-Over Factors.— amoment distribution meth- 
od, can partitioned yield the stiffness factors andcarry-over factors 
for the ends the member, The partitioning follows: 


rk 
| 
kips 
| 
at 
| 
bos 
les 
| 


whicheach and term canbe further partitioned, for example, follows: 


Eqs. 108 and 109, the terms designate moments and the terms designate 
forces the ends the member, The subscripts are read follows: 


represents the moment about the axis caused unit rotation 


about the axis with noother rotationor translationoccurring the ends 
the 
brief, the elements Eab represent stiffness factors and products 
stiffness factors and carry-over factors defined moment distribution, 
Reciprocal Relationships.—Much symmetry exists the construction 
because Maxwell’s reciprocal relationships for elastic members, For 
example, 


a 


a 
and 
a 


this instance, Maxwell’s reciprocal relationships are for reactions caused 
unit displacements, Note that the reciprocal reactions must correspond 
with the unit displacements accordance with the definition work. For this 
reason, 
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Sa)t # Sa)t” (117) 


NUMERICAL EXAMPLE 


Fig. shows structure consisting three prismatic members, each 
which parallel rectangular axis, The structure fixed ends and 
and subjected two concentrated loads. The assumed distributions mo- 
ments and forces are for structure cantilevered from end The structure 


TABLE 


Matrix 


-116 


and loads are the same those considered Figs. and previous 
work2 the writer and Michalos. 

The matrix analysis the structure summarized Tables and 
Each matrix inthese tables referred acorresponding equation inthe text, 
The matrices Table lare for determining the correction forces and moments 
any section the structure when the originof the and axes are end 
The magnitudes lg, and are the same those previously 
The signs however are different for several elements these ma- 
trices, because right hand rule axes used herein whereas left hand 


and 
r =10 0 1,50 Peo =10 -,833 62 
22.2 143, -412, 


rule was previously used the Table the correction moments 
and forces end are obtained, 

The matrices Table are referred the centroidal axes the structure, 
this table, the correction moments and forces the centroidal axes are ob- 
tained, 


TABLE 
Matrix Equation 
-10 -20 106,107 
-10 


Matrices and given Table for obtaining the stiffness ma- 
trix, Egp, the structure, The remaining elements can obtained 


means and Eqs. 104 through 107. 


CONCLUDING 
The present development has been extended tothe matrix analysis gen- 


eralized structure which lies space, and which continuous over several 
rigid flexible supports. 
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ANALYSIS CONTINUOUS TRUSSES CARRY-OVER 


Closure Jan Tuma 


JAN TUMA, writer appreciates the discussion Mr. 
Petcu, who pointed out the rapid convergence the carry-over procedure and 
analyzed the conditions The writer disagrees with Mr. Petcu’s 
statement, “that the carry-over procedure does not allow one obtain solu- 
tion with free terms appearing literally and that each load conditions must 
investigated separately.” The carry-over procedure may performed nu- 
merical (Part and algebraic form (Part 6). 

The algebraic form perfectly general terms algebraic carry-over 
factors and starting values. The general three and four 
recordedin Tables and The solution five and six span trusses re- 
corded Tables and 19, 

The special method2 suggested Mr. Petcu oneof many methods devised 
for the solution three-moment Its limitation seems the re- 
quirement equal spans. 

The writer indebted Oden, ASCE, who prepared the com- 
parative study four span truss the moment distribution method and the 
carry-over moment method (Tables and 21). Oden’s examples show the 
advantages and the simplicity the carry-over moment method, 


December 1959, Jan Tuma (Proc, Paper 2314) 
Prof, Civ. Engrg., Oklahoma State Univ., Stillwater, 
Cercetari Mechanica Aplicata,” Valeriu Petcu, Bucharest, No. 
October, 1959, pp. 855-865, 
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TABLE 20,—FOUR SPAN TRUSS, MOMENT DISTRIBUTION 


Analysis of a Four Span 
Continuous Truss by the 
Moment Distribution 


Method. 


(1) Angular Functions (Table 1): 


= 23.575 723 ° 19. 089 


(2) Moment Distribution Functions: 
(a) Stiffness Factors: 


1 
(>) Carry-Over Factors: 


. 978 
355 


Fixed End Moments: 
+ + 24.10 
k' 


FM' 


Distribution Factors: 


$ 
48' 48' 48' 48' 
* Kg, * Koy * | 
(3) Distribution Procedure: 
D +10 +10 
+ 3 it 
D + .05 |+ 06 | 
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TABLE SPAN TRUSS, CARRY-OVER MOMENT 


Analysis of a Four Span 
Continuous Truss by the 
Carry-Over Moment 


Procedure. 


(1) Angular Functions (Table 1): (>) Starting Moments: 


G Gin * Go3 G3 G5, = 0,355 


10 


710 " 732 23.575 19, 089 


Moment Carry-Over Functions: 


(a) Carry-Over Factors: 
CF, * CF, * OF, * Fig + Fig * 1, 956 


(3) Carry-Over Procedure 
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CONTINUOUS GIRDER BRIDGE 
WITH VARIABLE MOMENT 


Closure Sabri Sami 


SABRI ASCE.—Those who discussed this paper have presented 
valuable graphical and analytical alternate procedures. The discussions have 
added the value the paper and the kind comments are appreciated. 

The information presented Mr. Hayashi, ASCE, using the stiffness 
matrix analysis very interesting; actually the two methods are based the 
same principle. The main object the writer present simple method 
that any civil engineer can analyze easily and rapidly continuous bridge 
with variable moment inertia. 

The writerfails with Mr. Khachaturianthat the results the study, 
the form presented the writer, are academic interest only are ap- 
plicable very limited way. The writer believes that the paper presents 
practical and simple procedure for the analysis these bridges that the 
designer will able study structure with different relative dimensions 
using the simple formulas presented and thus choose the most economical di- 
mensions, The assumptions used the writer are clear and need not re- 
peated again; I-d (moment inertia, depth) curves can drawn and tenta- 
tive design canbe made for non-rectangular sections using the formulas the 
redundants presented the paper. 

The expressing the elastic equationof the redundant reaction the 
two-span bridge similar the approach Mr. Khachaturian,” and virtual 
work was used both studies. should noted that was Mr. Khachaturian’s 
encouragement that lead the attempt analyze the three-span bridge using 
the same approach, The writer expanded the method more than three spans 
and different rigid-frame bridges later 

Mr. Eremin, ASCE, has presented graphical method moment 
The writer grateful Mr. Eremin for renewing interest 
this interesting method, Fixed points can found easily using the analyt- 
ical results presented the paper. 

For example, the three-span bridge, case Fig. dividing Eq. 10b 
Eq. 10a one gets 


January 1960, Sabri Sami Paper 2346), 
Structural Engrg., West Virginia Morgantown, West Virginia. 
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This result can checked also from case (d) dividing Eq. 13b Eq. 13a 
follows: 


Therefore, two fixed points (because symmetry) can located Span BC. 
the four-span bridge, case (a), Fig. 13, dividing Eq. 15b Eq. 15a, yields 
The same results that can found dividing Eq. 21b Eq. case (d), 
Fig. 16. The minus sign has significance values are always nega- 
tive and K4-values are always positive. Therefore, fixed point can lo- 
cated each spans and CD. 
Dividing Eq. 15c Eq. 15b, case (a), Fig. Eq. 21c Eq. 21b, case 
(d), Fig. 16, the following obtained: 


(65) 
Mc Z4 


The same result obtained dividing Eq. 19c Eq. 19b case (b), Fig. 
dividing Eq. 22c Eq. 22b, case (e), Fig. 17. This locates another fixed 
point each Spans and CD. 

set fixed point used case the moment applied from the left 
from the right. The fixed points, one each end span and two each inter- 
mediate span, are independent the loading and are fixed location for any 
given bridge with different relative proportions. points can used, 
therefore, reduce the analytical work encourage graphical method, 
They indicate also the load position for maximum stresses, 

comparison between different straight parabolic haunches can stud- 
ied through the fixed-point This being done the writer and his 
students, The choice graphical semi-graphic solution depends the 
engineer’s general attitude towards graphical solutions, 

Mr. Leontovich has shown another approach the problem using method 
that based the concept elastic There are many approaches 
the problem published the United States and abroad; the writer’s method 
one them. The writer his method, Mr. Leontovich stated, 
“ready-to use” solution for determining the redundants continuous gir- 
ders parabolic shape. presented ina manner that can used easily 
and rapidly with the use ordinary office calculator. The aware 
Leontovich’s work and his valuable and important general method, 
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SECONDARY STRESSES PARALLEL WIRE SUSPENSION 


Discussion Janiszewski and Blair Birdsall 


and BLAIR BIRDSALL,? writers wish 
compliment the author the thorough-going attack has made subject 
that has had too little attention, After tests made 1922 and 1923 
preparation for the Philadelphia-Camden Bridge that producedsome empirical 
results, designers have been content assume that the subject secondary 
stresses were sufficiently inconsequential require more attention than 
little judgment setting the factor safety for bridge cable. Thus, thor- 
ough analytical study very muchin order meansof either confirming 
revising these empirical conclusions, 

Unfortunately, the press every-day events have denied the writers the 
pleasure and satisfaction full review and detailed mathematical critique 
this paper. However, they hope that the few remarks they have been able 
piece together will stimulate the author close the discussion correlating 
his mathematical speculations with the results known tests and with exper- 

For the benefit those who are more mathematically sophisticated than 
the writers, they offer the following description their struggle overcome 
apparent paradox the text: 

dealing with beam problems, one accustomed associating greater 
stiffness with greater unit stresses for given deflection. Thus, the writers 
were puzzled when, studying Eqs. and they found that, according 
these expressions, unit secondary stress (Eq. increases stiffness de- 
clarify, note that although defined only “flexural stiffness 
one panel cable,” applied mathematically being equal Be- 
cause (in which the metallic area the cable), this explains the 
presence the expression for becomes smaller, increases. 
secondary stress increases, Therefore, stated previously, stiffness de- 
creases, unit bending stress increases, 

After first falling into the “trap” there was 
the author’s reasoning and thereby following dead-end trail, the writers 
returned take closer look the reasoning, and last found the solution 
the apparent paradox: 

For given lateral displacements certain points abeam with respect 
others, unit stress does increase with stiffness. However, are given re- 


July 1960, Thomas Alan Wyatt (Proc, Paper 2542), 

Senior Design Bridge Div., John Roebling’s Sons Div., The Colo. Fuel and 
Iron Corp. 

Engr., Bridge Div., John Roebling’s Sons The Colo, Fuel and Iron Corp. 
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lative rotations certain cross sections abeam with respect others, unit 
stress independent stiffness. The rotations simply force changes length 
the fibers. Normally, therefore, such case, change would not af- 
fect the secondary stress. this particular case, however, the rotation it- 
self function the stiffness. stiffness decreases, the relative rotation, 
and thereby, the secondary stress, increases. Placing this thought the prac- 
tical context the subject problem, see the perfectly logical inter-relation: 
let consider the end panel the cable adjacent tower the 
cable relatively flexible, the rotation caused the tower saddle changes 
load absorbed large extent the first panel, and the forces tending 
rotate the cable band the endof panel are relatively small, Thus, 
the cable band rotates only and the relative rotationof the ends 
the first panel cable almost equal the entire rotation the cable 
cross-section the saddle. 

the other hand, the cable relatively stiff, the cable will force the 
first cable band rotate, thereby transmitting some the rotation into the 
next panel, 

considering the models presented through has occurred 
the writers that another method presentation might have made the anal- 
ogy clearer many readers, 

The author makes the simplifying assumptions that, before end rotation 
applied, the cable, any panel point straight through the cable band, that the 
suspender tension does not exist, and that cable tension does not change the 
band, then proceeds express the moment existing any panel point af- 

Another way illustrate this phenomenon consider the cable polygon 
is, the cable angle and tension changing each panel point, the pres- 
ence the suspender tension, The common assumption that the cable repre- 
sents string polygon implies that the three forces each point, two cable 
tensions and suspender tension, meet point locatedin the line the sus- 
pender force. One could then reason that this string polygon disturbed 
the secondary influences which the author studying. This disturbance has 
the affect moving the point intersection the cable tension slightly to- 
ward the tower toward mid-span distance (for eccentricity). The ver- 
tical reaction produced the cable tensions this new point intersection 
would still equal the suspender pull (say The moment would then 
expressed 

The writers not doubt that the author’s approach lends itself better 
mathematical manipulation, but the foregging may help clarify the physical 
concept, 

Experience this subject could probably summarized saying that 
conditions have been observed that have given designers cause for concern. 
the same time, many designers recognize the existence secondary stresses 
this type and consciously attempt minimize rotation, One reason for the 
lack observable affect may the capacity cold-drawn hot galvanized 
bridge wire absorb great over-stress with little permanent elongation and 
with resultant improvement elastic properties. 

few examples from practice may interest: 

Persons concerned withthe design and the Forth Road Bridge 
Scotland noted that one the saddles relatively large rotation 
between free-hanging cable and dead-load conditions, conventional proce- 
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DISCUSSION 


dures were followed, minimize this rotation, means were found accom- 
plish the necessary motion partly rotation and partly 

canvassing experience prior reaching this decision, the following points 
interest were developed: 

During construction the San Francisco-Oakland Bay Bridge, one the 
saddles experienced rotation equal that which gave rise the discussion 
the Forth Road Bridge. noticeable ill-effects have been observed. 

similar saddle the Mackinac Bridge, the designers purposely made 
the member supporting the saddle long possible, minimize the rota- 
tion, 

The writers know instance where broken wires have been found that 
could attributed the subject Careful inspections have been 
made many existing bridges without uncovering any problems this type. 
One the most complete such surveys was that recently conducted the 
Brooklyn Bridge, preparation for its new life modern vehicular bridge. 
Over the years, wiresfrom various locations the Brooklyn Bridge 
cables have been tested and found have properties good, not better, 
than their original state. 

However, indicated previously, specific design criteria are always bet- 
ter than rules-of-thumb pure judgment and the author encouraged carry 
his thinking through that point. 

source confirmation correction his reasoning, the author 
referred areport dealing with the Delaware River Bridge whichthere will 
befound data having adirect bearing onthe this paper.4 Correlation 
these results with the mathematics the paper would great value. The 
author may gratified note the experimental verification contained therein 
the series cable “kinks” shown the cable bands Figs. and will 
also noted that the experiments developed friction coefficient between 
15% and 20% between the wires wrapped 

Under the heading “Slip Through Hanger Bands” the author co- 
efficient friction 60% has been suggested for the friction betweenthe wires 
within cable band, may that this value the result erroneous in- 
terpretation test results. Tests pertinent this point are also included 
the above-mentioned Here seen that the force required move 
cable band along parallel wire cable was approximately 60% the total 
bolt tension, Several other tests other projects have confirmed this value. 
should noted, however, that two friction surfaces are involved, There- 
fore, the actual coefficient friction one-half this ratio, approximately 
30%. These figures are artificially high with respect other pure friction 
tests, The difference generally attributed the additional resistance 
sliding caused the bulging the cable adjacent the band, Thus, there 
some reason question the use even the 30% value appropriate co- 
efficient friction for this study. 

The same report® shows the force required pullindividual wires out 
sample cable clamped cable band, this case, the value was aver- 
age corresponding tothe 60% value stated previously. Thus, the cor- 
responding true coefficient friction would approximately This leads 


Report the Board Engineers the Delaware River Bridge Joint Com- 
mission,” pp, 124 and 125, 
pp. 126 and 127, 
127. 
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one wonder the range 15% 20% should not considered appropriate 
for the wires under the band well for wires under wrapping. 

The writers hope that, for his closing discussion, the author will find time 
use, least comment on, the suggestions contained herein, correlate 
his reasoning with the available test results, and conclude with some spe- 
cific design criteria guide the designer. 
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DYNAMIC RESPONSE ELASTO-PLASTIC FRAMES* 


Discussion John Blume 


JOHN ASCE.—This paper notable addition the still 
meager fundof information the inelastic response structures idealized 
dynamic systems actual earthquake records. With continued efforts this 
direction sufficient data will some day hand isolate and assign re- 
liable numbers the many parameters the problem. 

may postulated that the two basic items attenuation for the author’s 
single mass elasto-plastic system are damping and work done inealstic re- 
However, each these factors has subdivisions and there doubt 
interaction between the sub-divisions and the basic items, well with res- 
onant magnification. Changes the displacement ratios (Figs. 11, inclu- 
sive) can considered indicating the values the broad categories 
damping and inelastic response within the limits the author has noted, The iso- 
lation the two attenuating factors seems important for bothresearch and in- 
elastic energy design 

Mr. Penzien does not present data for the cases infinite yield point with 
damping 0.05, 0.10, and Such would help isolate 
the damping attenuation any period the elasto-plastic attenuation 
obtained deducting that due damping from the total. This can 
done directly withthe data provided forthe zero damping condition and the re- 
sults are most interesting, The author’s data indicate that the effectiveness 
the inelastic response varies (when with not only but with period, 
The reasons for this doubt involve the proximity the quasi-resonant 
magnification peak the short period range the Centro 1940 earthquake 
spectrum, 

factor sometimes overlooked this fieldof effort that designer may 
consciously unknowingly change stiffness with strength can al- 
change the yield value without changing stiffness changes, 
does period, assuming the mass essentially constant. Should proposedstruc- 
ture initially designed have 1,2 sec fundamental period initial 
yield value 0,12 both strengthened and stiffened (in the later design 
stage) that became 1.0 sec and yield 0.15 the response attenuation may 
essentially unchanged even The reverse may true, however, for 
lower yield values. 

Multistory buildings, especially the skeleton and skin type without walls, 
may well have lower-story yield values 0.10 and less. They probably have 
damping values, distinguished from inelastic energy absorption, the or- 


July 1960, Penzien (Proc, Paper 2545). 
Pres., John Blume Assoc., Engrs., San Francisco, Calif, 
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der critical, and certainly not over 10%, Further inelastic studies 
this range would most interesting. 

study the author’s figures indicates several matters and raises ques- 
tions for future research including the following: 


any period peak the spectrum, the relative at- 
tenuation values damping vary greatly with the 

Damping value the shorter periods relatively more effective for 
greater than for low yield points. 

The relative difference between 5%, even 15%damping any 
period for all except very low yield values minor; that is, general 10% 
damping offers minor attenuation over that for damping. 

damping provides, especially for large elastic response, 
considerable reduction over damping for the entire period range. 

For given period and damping value, yielding inelastic condition 
tends attenuate response otherwise (elastically) greater magnitude than 
the yield value, except for very low yield values, below which the response in- 


These and other indications are general agreement with other test data 
and proposed approaches inelastic design. 

The “optimum” yield values which reference made Conclusions 
and should approached with caution, Mr. Penzien notes, until other 
earthquake records are tried. The N-S component the same earthquake, 
for example, generally more severe than the E-W component utilized this 
work, Other earthquakes variable epicentral distances, intensity, and soil 
conditions must also considered, noted that for damping values 
other than the optimum not true optimum but essentially constant for 
greater than the 0.10 gto 0.20 range. This indicates elastic re- 
sponse for this earthquake since the viscous damping assumed was sufficient 
provide dynamic stability without inelastic excursions. interesting ques- 
tion how would damping perform for the stronger N-S component 
Centro 1940? 

The relatively small range from code design unit stresses yield unit 
stresses combined with code base shears and the high accelerations elastic 
earthquake spectra plainly indicate that inelastic response basic part 
the earthquake resistance most buildings.3 Damping alone cannot dothe job. 
More work inelastic response, (elastic-plastic and other types) with many 
earthquake records, variable stiffnesses, and with multi-mass systemsas well 
single mass systems desirable. The profession also needs data 
from more static testing members and frames, not just yield but with ade- 
quate readings all the way failure, and with reverse loadings and repeated 
cycles having various degrees inelastic deflection. take “our heads 
out the sand” see the economic need for emergency inelastic 
earthquake response shall make progress, doubt much more than the 
fabled ostrich, 


TESTS 120-FT SPAN PRESTRESSED CONCRETE 


EREMIN,? this paper Base and Rowe have 
shown interesting design and test measurements con- 
crete beam, that are similar those used the overpass the Chiswick 
Rotary Junction the Great West Road, England, 

The authors have proved that the assumptions used design the bridge 
beams were close agreement with the resulting test measurements, How- 
ever, the study the physical and elastic properties the test beam, 
interesting note the statement made the authors’ the result- 
ing tests measurements (under the heading, “Analysis and Discussion Re- 
sults: Test Results”): “after the first twocables had been stressed, the move- 
ment the joints continued until five six cables had been 

Appearance cracks the joints and their closure, only after five six 
cables the beam were stressed, shows that the details the joints require 
further improvement, 

Furthermore, Messrs. Rowe stated that the east joint was opening 
lead approximately kips per jack whereas the midspan the crack 
actually commenced slightly under 134 kips. This fact also indicates that 
the careful consideration for further improvement design and method 
casting the joints between precast sections needed, 


September 1960, Base and Rowe Paper 2608). 
Cons, Structural Engr., Berkeley, 
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FREEZING AND THAWING EFFECTS PRESTRESSED 


Discussion Walter Kunze 


WALTER KUNZE, ASCE.—One the conclusions stated inthis pa- 
per durability improved Although there are ample 
data substantiate this conclusion, the data presented Table not nec- 
essarily lead this The average durability factor the 3,000 psi 
plain concrete shown for six samples, Omitting specimen D-2 (that 
53.7 points lower than any the other five) the average becomes that 
higher than any the average values for 5,000 psiconcrete and thus was 
more durable than the prestressed concrete higher strength, This might 
due air entrainment the mix whereas other mixes may have been non 
The paper indicates that air-entraining admixture was used 
but from Table only mix IID appears substantiate this. From this stand- 
point, the authors seemto emphasize the need for air entrainment producing 
durable concrete whether not prestressed, 

The paper clearly shows that stress level 2,000 psi compression, 
prestressed concrete more durable than unstressed concrete, Prestressed 
concrete commonly usedin flexural members, Although the maximum pre- 
stress member isoften 2,000 psi higher, the average far Fur- 
thermore, under design load, much the concrete stressed very low 
level and some This again points out the need for durable con- 
crete, whether not the concrete For most concretes used 
resistance tofreezing and thawing, along withlow water-cement ratio and ade- 
quate curing. 


October 1960, Gutzwiller and Musleh (Proc, Paper 2633), 
Mgr., Structural and Railways Bur. Portland Cement Assoc. 
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PRINCIPLE VIRTUAL WORK STRUCTURAL 


Discussion Ray Clough 


RAY CLOUGH, has been presented interesting 
paper demonstrating the general applicability the principle virtual work 
problems instructural mechanics, However, his not point 
out the dual nature the principle, that, the opinion the writer, the key 
its most effective utilization. The examples presented Mr. DiMaggio rep- 
resent both aspects the duality, but the distinction between the two forms 
which the principle may applied not made clear, Although the dual nature 
the virtual-work principle has been described before, brief restatement 
may interest. 

The basic principle virtual work has been expressed mathematically 
Eq. words this principle may stated follows: 


“If any structure equilibrium under the action set exter- 
nal forces and subjected displacement which compatible with 
the constraints the structure, the work done the external forces 
during this displacement equal the work done the internal forces.” 


The word “virtual” the this work principle denotes the fact that the 
mechanism causing the displacements not associated with the forces acting 
the structure, 

The validity the principle may demonstrated, Mr. DiMaggio has 
shown, starting with the equations equilibrium. Conversely, the work 
principle could accepted fundamental and the equations equilibrium 
derived therefrom. Obviously, makes little difference which point view 
adopted. The important point note that inherent the application 
the principle virtual work are two requirements: The system must 
equilibrium and the displacements must compatible with prescribed con- 
dition contraint. The forces and displacements need not (and normally 
are not) related effects, however. 

The dual nature the virtual-work principle results from the two ways 
which may applied, the system subjected set forces equi- 
librium, the expression the equality internal and external virtual work 
will impose the requirement compatibility. other hand, the sys- 
tem subjected compatible deformation pattern, the equality internal 
and external virtual work will express the requirements equilibrium, Fol- 
lowing procedures, the principle virtual work may eval- 
uate either displacements forces structural system. 


November 1960, Frank Maggio (Proc, Paper 2643), 

Civ. Engrg., Univ. California, Berkeley, 

ing, Vol. XXVI, October, 1954, 
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Considering the first use, the basic problem isthe evaluationof certain de- 
sired displacement quantities resulting from given condition deformation 
that exists inthe structure. determine displacement component that 
compatible with these deformations, artificial (or virtual) force system 
applied. The only requirement imposed these virtual forces that they 
must equilibrium; they need bear relation, whatsoever, any actual 
system forces reactions acting the structure. Making use this ap- 
proach, the force quantities Eq. wouldrepresent virtual force 
system equilibrium, while would the actual deformations and com- 
patible displacements. The mechanism causing these deformations unim- 
portant. Mr. DiMaggio has pointed out, they may elastic deformations 
due actual loads, they might visco-elastic effects, they might caused 
temperature changes, and forth. However, this use the principle, 
should clearly understood that these deformations are definite, actual ef- 
fects, and they should not termed The force system the virtual 
system this 

the second use the virtual force principle, the objective the evalua- 
tion certain specified forces (usually external) that are equilibrium with 
given (usually internal) force pattern, determine the desired force com- 
ponents, appropriate artificial (or virtual) system displacements in- 
troduced, The requirement imposed these virtual deformations that they 
must compatible with the essential constraints the system. However, they 
need not satisfy all the actual support internal continuity conditions, 
fact, displacement must introduced deliberately any point where aforce 
determined; and the desired quantity internal force re- 
action, the displacement will violate actual constraint (as the last ex- 
ample). Making use this approach, the deformations n(x) Eq. would 
represent virtual displacement, while the forces (W, would ac- 
tual equilibrium force These forces might due any cause: load- 
ing, thermal changes, displacements supports, and forth; but they must 
actual forces known amplitude. 

his examples, Mr. DiMaggio has illustrated basic uses the vir- 
tual work principle. also has correctly demonstrated the freedom choice 
that inherent the selection virtual force displacement patterns, that 
is, the fact that the virtual force system independent the actual forces and 
reactions inthe structure, and that displacement system not con- 
trolled continuity requirements the actual However, the speci- 
fic requirements the principle virtual work which are summarized pre- 
viously have not been clearly stated Mr. DiMaggio. 

Furthermore, the statement under the heading “Restricted Deviation for 
Beam”: 


“The function (x) does not necessarily satisfy any the constraints 
deflection the structure acted W.” 


true only whenthe principle virtual displacements applied, but not true 
the principle virtual For example, any Mr. DiMaggio’s ex- 
amples deflection computations (that are applications the principle 
virtual forces), essential that n(x) the actual displacements the 
structure, satisfying all boundary conditions and constraints the portions 
the structure which deflections are computed, 
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One other word caution offered, Eq. 57, the curvature term due 
shear forces valid only for prismatic members; otherwise, the 
variations and with would lead additional terms the expression. 
the general case straight beam with arbitrarily varying section, the 
correct expression for deflection flexural element (including the effect 
shear deformation) given Eq. but not Eq. 59. 
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PERIODS FRAMED BUILDINGS FOR EARTHQUAKE 


Discussion Merit White 


shows how easily one can take account various kinds flexibility (shearing 
and bending the structure, and rotation and sidesway its foundation) 
computing the lower natural periods vibration building. However, two 
limitations its practical usefulness may mentioned: (1) the structures 
that have the characteristics assumed the derivation Eqs. 28a and are 
rare, and (2) the parameters that appear these equations are easier de- 
fine than may added that the writer these equations seem 
rather complicated use. 

The structure that envisaged the authors has the following character- 
istics: 


(a) The weight per unit height per story constant throughout the struc- 
ture except the roof that heavier than the floors. Thus, the structure 

(b) The shear stiffness any level directly proportional the total 
weight above that level, except the first has just half the stiffness 
the second story; this Because the contributionof the struc- 
tural frameto the stiffness normal very small compared that 
its walls and partitions unlikely that the shear stiffness structure 
uniform section will increase factor order between and infinity 
descending from the roof the second story would have re- 
main proportional the total weight above, and implied the values 
assume constant shear stiffness. 

(c) The flexural stiffness increases amount proportional distance 
below the roof. The validity this assumption depends onthe extent which 
walls and partitions not participate resisting bending. 

(d) The soil supporting the structure allows both rocking and lateral dis- 
placement proportional tobase moment and base shear, respectively. Finally, 
the natural period vibration corresponding each these types distor- 
tion determined and the period vibration due all them combined 
found squaring, adding, and taking the square root the sum. 


The authors verify their results computing the lower mode periods 
the Alexander Building San Francisco utilizing the structural parameters 
(weight, and forth) published Because the com- 
puted periods agree fairly well with those observed the authors conclude that 
their relations are reliable and can used for accurately predicting periods 


December 1960, Mario Salvadori and Ewald Heer Paper 2680), 
Head, Civ, Engrg., Univ. Massachusetts, Amherst, Mass, 
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This true far the relation between structural param- 
eters andperiods However, the accuracy the result limited 
the accuracy with which the parameters are known, The stiffness param- 
eters quoted were obtained not measurement direct com- 
putation but indirectly, trial and error computations periods and mode 
shapes that were repeated with different parameters until the predictions agreed 
with observations. necessary find periods from examination 
structural plans, the accuracy will worse, although probably better than that 
Eq. 

far the actual computation periods from given set parameters 
concerned the method given herein probably one the simplest and 
the sametime There restrictions mass and stiff- 
ness distributions. 

First, reasonable deflection curve chosen, say concave inward para- 
bola shear the principal source concave outward parabola 
bending principal importance, astraight lineif both equally 
significant base rotation appears subsequently, does not 
make very much difference which one selected. Next, the distribution 
maximum accelerations tabulated plotted; simply the displacement 
each point multiplied square the yet unknown angular fre- 
quency The maximum acceleration any section multiplied the corres- 
ponding mass per unit height structure gives the inertia loading per unit 
height. From the inertia loading the shear (V) and the bending moment (M) 
along the lengthof the building graphical integration from the roof 


where the shear stiffness defined the authors and the bending 
stiffness, The totalelastic energy the structureis the area beneath the elas- 
tic energy distribution curve. this must added the elastic energy stored 
the foundation that equals lateral displacement neglected, 
usually permissible. The elastic energy the structure-foundation sys- 
tem will contain the unknown common factor 

The maximum kinetic energy next found, The assumed displacement curve 
multiplied givesthe velocity distribution. This squared and multiplied 
m/2 gives the kinetic energy distribution along the building mass per 
unit height and need not uniform). The area beneath the kinetic energy dis- 
tribution curve the total kinetic energy the structure, This will contain 
the unknown When this system vibrates one its natural modes there 
cyclical transformation energy from kinetic elastic kinetic, and 
forth. that the maximum elastic energy equal the maximum kinetic 
energy. Equating the two energies just computed gives expression contain- 
ing the single unknown p2, that determined thereby. 

This method becomes much simpler and still reasonably accurate for 
many structures the mass and the stiffnesses are allassumed constant 
and equal the average values the actual quantities. The following sym- 
bols are introduced addition modifying those defined the authors: 


maximum displacement the roof (in feet) 


inertia force per unit height (in kips per foot) 


average cross-sectional moment inertia for flexure (in ft4) 
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average shearing stiffness coefficient (in kip-feet per foot) 


bending moment any section (foot-kips) 


maximum velocity any section (feet per second) 


Assuming the extreme displacement the fundamental mode trian- 
gular, zero ground and roof level, one gets the equation motion 
any element 


and the maximum velocity 


The shear and bending moment any level due the inertia loading are 
and 
The moment the foundation found letting Eq. 34, 


the lateral displacement the foundation neglected, the total elastic en- 
ergy the system 


(36) 


The total kinetic energy 


Equating the two energies and solving for gives 


This corresponds Eqs. 28a and 29, 
For deflection assumed concave outward parabola the expression 
for the period becomes: 


Then the maximum inertia loading 
4 
(39) 


For quarter sine curve, convex outward, one obtains 


When these expressions are applied the north-south fundamental mode 
the Alexander Building, using the following parameters: 2,2, 196, 
ively, 1.29 sec, 1.26 sec, and 1.28 sec. 

The observed periods for this mode ranged from 1.27 sec. 

The writer considers this method period determination have several 
advantages over the one described the authors: (1) the assumption con- 
stant more realistic for the great majority uniform buildings 
than the assumption that shearing stiffness proportional tothe total mass 
above any section, (2) the expressions that have been obtained for the period 
uniform structure are simpler use than one that requires 
reference tables for values coefficients, (3) the method can easily 
generalizedand applied structures that nonuniform require spe- 
cial treatment. 
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Journal the Structural Division 


Proceedings the American Society Civil Engineers 


“Load Distribution Composite Girder-Slab Systems” Hondros and 
Marsh, November 1960, Paper 2645. 


The caption for Fig. should follows: 

FIG, AFTER CENTRAL LONGITUDINAL CUT 
81. The caption for Fig. should follows: 

FIG. FORMWORK AND REINFORCEMENT 


Eq. the terms the first set parentheses should fol- 


“Movements Cable Due Changes Loading” James Michalos and 
Charles Bornstiel, December 1960, Paper 2674, 


pp. and 35. Table should have its two halfs interchanged that the 
half headed TABLE page and the half headed TABLE 5,—CON- 
TINUED appears page 35, 


ore 
wht, 
7 
Je q 
: 
ta 
= 


= 
: 
j 
: 
2 
: 


PROCEEDINGS PAPERS 
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